EXAFS study of lead-free relaxor ferroelectric BaTii^Zr^Os at the Zr K-edge 
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Extended X-ray absorption fine structure (EXAFS) experiments at the Zr K-edge were carried 
out on perovskite relaxor ferroelectrics BaTii-zZr^Oa (BTZ) (x = 0.25, 0.30, 0.35), and on BaZrCb 
for comparison. Structural information up to 4.5 A around the Zr atoms is obtained, revealing 
that the local structure differs notably from the average Pm3m cubic structure deduced from X-ray 
diffraction. In particular, our results show that the distance between Zr atoms and their first oxygen 
neighbors is independent of the Zr substitution rate x and equal to that measured in BaZrOa, while 
the X-ray cubic cell parameter increases linearly with x. Furthermore, we show that the Zr atoms 
tend to segregate in Zr-rich regions. We propose that the relaxor behavior in BTZ is linked to 
random elastic fields generated by this particular chemical arrangement, rather than to random 
electric fields as is the case in most relaxors. 

PACS numbers: 61.10.Ht, 77.80.-e, 77.22. Gm, 77.84.Dy 
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I. INTRODUCTION 

The perovskite-type barium zirconate titanate, 
BaTii-zZr^Oa (BTZ), has attracted considerable at- 
tention as a possible lead-free ferroelectric material 
to replace the current industry standard lead-based 
ferroelectrica 1 ! 2 ! 3 ! 4 ! 5 ! 6 . The particularity of BTZ is the 
continuous change of its properties from a ferroelectric 
behavior at low Zr substitution rates to a relaxor 
ferroelectric (relaxor) behavior at higher substitution 
rates (0.25 < x < 0.5)£. The low substitution regime 
has been mainly studied for its interest in lead-free 
ferroelectric memories^. However, due to the recent 
report of outstanding electro-mechanical properties in 
relaxor-based ferroelectric solid solutions^, it is rather 
the relaxor regime of BTZ which currently attracts a 
large research effort. 

Relaxors are characterized by a broad and frequency- 
dependent dielectric anomaly as a function of tempera- 
ture, instead of a sharp and frequency-independent di- 
vergence as in classical ferroelectrics. In a relaxor such 
as PbMg 1 / 3 Nb 2 /303 (PMN), the anomaly surprisingly 
occurs in a cubic average structure which remains cen- 
trosymmetric (non polar) at all temperatures. It is 
generally admitted that the peculiar properties of re- 
laxors are related to the presence of nano-scaled polar 
regionsiSiii* 1 ^, due to different cation shifts in different 
parts of the structure. 

Two necessary ingredients are often cited for a re- 
laxor ferroelectric: (i) the presence of Pb 2+ or Bi 3+ 
cations (showing large displacements due to their lone- 
pair) and/or (ii) a heterovalent cationic disorder (gen- 
erating random local electric fields 1 ^ that break long- 
range polar correlations). They are both present in the 



extensively studied model relaxor PMN (e. g. Refs. 
[i^lTMl^lT7lli^lT9ll2^ . On the other hand, BTZ relax- 
ors present a homovalent Zr 4 +/Ti 4+ substitution, which 
does not give rise to such random electric fields. There- 
fore, another mechanism has to be considered in order to 
account for the break of long-range correlated displace- 
ments in BTZ. It has been proposed that random elastic 
fields, induced by the difference in size of Zr 4+ and Ti 4+ 
cations, play an important role2i. Yet it is not known 
how, and on which scale such random elastic fields build 
up. The BaTii-^Zr-jOs solid solution thus opens an in- 
teresting route towards the comprehension of the relaxor 
behavior, and raises much interesliSiS 2 * 2 ^^ 2 ^ 2 ^. Further- 
more, it is intriguing that PbZi'i-^TLjOs (PZT) shows 
no relaxations 7 , while BTZ does. 

The average crystal structure of the BTZ relaxors 
(0.25 < x < 0.5) is cubic (space group Pm3m) at any 
temperatureS 4 . Local distortions away from the ideal cu- 
bic structure expected from the observation of the re- 
laxation behavior are evidenced by the observation of a 
first-order Raman scattering 2 ^ 3 ., which is forbidden in 
a perfect primitive cubic structure. However, the latter 
technique does not reveal the nature of these distortions. 
The aim of the present study is to determine the local 
structure and chemical order in BTZ relaxors. For this 
purpose, we used EXAFS (Extended X-ray Absorption 
Fine Structure) spectroscopy which is a structural and 
chemical local probe. This technique has been proven 
to be a powerful tool to analyze the local distortions in 
ferroelectric perovskites (e. g. Refs. I2cl29l30ll3ll) . On 
the other hand, only a few EXAFS studies have been 
reported for relaxors yet 32 : 33 ' 34 : 35 ' 36 . 

In the following, we report an EXAFS study of three 
BTZ relaxors (a: = 0.25, 0.30 and 0.35) at the Zr K-edge. 
The local structure parameters were determined up to 
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distances of 4.5 A from the Zr atoms, giving information 
on their local environment up to their fourth neighbors. 
Hence, the Zr/Ti repartition could be probed, as well as 
the distortions induced by the substitution of Ti atoms 
by Zr ones in BaTiOs. In order to have a proper reference 
for the understanding of the local Zr environment, we also 
studied BaZr03 under the same experimental conditions. 



II. EXPERIMENTAL METHODS 
A. Samples and experiments 

BaTii-^Zr^Oa (x=0.25, 0.30, 0.35 and 1) powders 
were synthesized by solid state reaction, starting from the 
appropriate amounts of BaC03, Ti02, and ZrC>2 powders 
and following the method described in Ref. 7. They were 
then characterized by X-ray diffraction at 300 K. All the 
powders were found to be single phases within the accu- 
racy of the experiment. The cubic cell parameters are 
4.054(1), 4.061(1), 4.072(1), and 4.192(1) A for x values 
of 0.25, 0.30, 0.35, and 1 respectively. 

EXAFS experiments at the Zr K-edge (17.998 keV) 
were carried out on the FAME beamline BM30B at the 
European Synchrotron Radiation Facility (ESRF). The 
X-ray absorption coefficient \i was measured in the trans- 
mission mode as a function of the incident photon energy 
E in the range [ 17.8-19.5 keV ], using a Si(220) single 
crystal monochromator. The samples consisted of pel- 
lets made of BaTii-^Zr^Os powders (with a grain size 
of about 1 (im) mixed with very low-absorbing boron 
nitride. The effective thickness e of each sample was 
chosen in such a way that the product /xe increases ap- 
proximately by 1 through the Zr K-edge. EXAFS data 
were collected at room temperature for all samples, as 
well as at 11 K for BaZr0 3 and at UK, 90 K and 150 K 
for BaTio.65Zro.35O3. 



B. EXAFS data treatment 

The scattering of the photoelectron by the neigh- 
bors of the central absorbing atom introduces oscilla- 
tions in the energy dependence of the absorption coef- 
ficient /J,(E) after the energy edge. The normalized ex- 
perimental EXAFS signal is given by x(fc) = [n(k) — 
fio(k)]/ Afio(k — 0). jLto is the smooth, atomic-like, ab- 
sorption background, A.fj,o(k = 0) is the absorption edge 
jump and k is the photoelectron wave number, given by 

k 



2m e (E — Eo)/h with m e the electron mass and 
Eq the edge energy. In the present work, x(^) was ex ~ 
tracted from the measured absorption coefficient (J,(E), 
by using the AUTOBK program^. E was taken at the 
first maximum of the derivative of n(E) at the Zr K-edge. 
Two representative examples of normalized fc 2 -weighted 
x(fc) spectra are shown in Fig. ^ 

The theoretical EXAFS signal x(k) is expressed as a 
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FIG. 1: fc 2 -weighted normalized EXAFS signals, for BaZr0 3 
and BaTio.65Zro.35O3 at 300 K. 



sum of contributions from different paths, each path i 
corresponding to a given scattering process of the photo- 
electron: 



X (k) = -s 2 J2 



2 ^ NjAj(k) 2 fcV 2 -2Rj/x 



(fc) 



kRi 



sin (2kR, + 25 c (k) +&(*;)) ; 



(1) 



where Ni is the degeneracy of path i, Ri its half-length 
and Ai(k) its effective scattering amplitude. The Debye- 
Waller (DW) factor a 2 is the standard deviation of the 
Ri distance distribution, assumed to be gaussian. The 
DW factor takes into account both the thermal disor- 
der and a possible small structural disorder. S c (k) and 
4>i(k) are phase shifts associated with the electron prop- 
agating into and out the potentials of the absorbing site 
and scattering sites respectively. The other parameters 
are the photoelectron mean-free path A(fc) and an over- 
all amplitude factor Sq 2 , close to 1, which accounts for 
many-electron effects in the excited central atom. Note 
that Eq. ^ includes both single (back-)scattering (SS) 
and multiple scattering (MS) processes^. For a SS path, 
Ni is simply the number of chemically identical atoms 
situated at a given distance Ri from the central atom. 

The analysis of EXAFS signals at the Zr K-edge was 
performed in i?-space after a Fourier transform (FT) 
of the fc 2 -weighted x(fc) in the fc-range [ 3.2-14.6 A -1 ], 
using a Hanning weight function. The calculated FTs 
were fitted to the experimental ones using the FEFFIT 
program^. The electronic parameters Ai{k), <f>i(k), S c (k) 
and X(k) were calculated for both SS and MS paths with 
the FEFF8.00 code^. The structural parameters Ri, a 2 , 
as well as Ni if unknown, were extracted from the fit as 
will be described in details in sections IIII Bl IIII CI and 
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IIII Dl As customary, an additional parameter AEq was 
introduced to account for the small difference between 
the experimental edge energy E and its calculated value 
using the FEFF code. The parameter Sq 2 was also re- 
fined, the experimental value usually differing from the 
theoretical one. 



III. RESULTS 
A. General trends 

Examples of the Fourier transforms of the experimental 
fc 2 -weighted x(k) spectra are shown in Fig. [21 F° r all 
studied samples and at all temperatures, the FTs can be 
readily separated into two contributions, characteristic 
of the perovskite structure. The first neighbors of a Zr 
atom are six oxygen atoms (denoted hereafter as 01), 
which form a regular octahedron in the perfect Pm3m 
perovskite structure . The first contribution to the FT, 
at R values lower than 2.5 A, only includes a SS process 
by one of the 01 atoms. The second contribution, in the 
i?-range [ 2.5-4.5 A], results both from SS processes where 
the photoelectron is backscattered by either the second 
(Ba), third (Zr or Ti), or fourth oxygen (02) neighbors of 
the Zr atom, and from several MS processes. Note that 
due to the /c-dependence of the phase shifts 4>i(k) and 
S c (k), the maxima of the FT modulus occur at distances 
different from the real ones. 

Surprisingly, the first peaks of the FTs are nearly iden- 
tical at 300 K in all relaxor samples and in BaZr03 (see 
Figs. |21 and The latter suggests very close first neigh- 
bors environment of the Zr atoms, whatever the Zr sub- 
stitution rate. On the other hand, the FTs of BaZr03 
and BTZ relaxors differ significantly for R values larger 
than 2.5 A. 

The shape of the FTs suggests an analysis in two steps. 
Fits of the measured FTs are first performed in the R- 
range related to the Zr first oxygen neighbors (Sec. IIII B|l . 
In a second step, the fitted i?-range is extended up to 4.5 
A, in order to take into account the further neighbors 
contributions (Sec. IIII CI and IIII Dfl . In BTZ relaxors, 
one has to take into account both the Zr/Ti substitution 
and the strong contribution of several MS paths in this R- 
range, which makes the EXAFS analysis complicated. In 
order to get reference parameters for the analysis of BTZ 
relaxors, the FT of BaZr0 3 is first fitted (Sec. ITTTT^I . 
From previous EXAFS studiea 41 i 42 i 43 , the local structure 
of BaZr03 can be considered as a perfect cubic perovskite 
structure 4 ^, identical to the average one deduced from X- 
ray diffraction (space group Pm3m)^i&. 



B. Analysis of the first neighbor contribution 

We already noted that the backscattering processes be- 
tween Zr atoms and their first oxygen neighbors (01) 
give rise to very similar contributions to the FTs of the 
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FIG. 2: Selected FTs of the k -weighted EXAFS signals: 
BaZr03 and BaTio.65Zro.35O3 at 300 K (top), and thermal 
evolution in BaTio.65Zro.35O3 (bottom). The imaginary part 
and the modulus of the FTs are plotted as dashed and solid 
lines respectively. 



EXAFS signals below 2.5 A, for all samples and at all 
temperatures (Figs. [2land[!3)- From the known structure 
of BaZrOa, it is then reasonable to assume, as a starting 
point for the analysis of BTZ relaxors, that the six 01 
atoms are located at the same distance from the absorb- 
ing Zr atom. Within this hypothesis, only one SS path 
(denoted hereafter as path 1) contributes in Eq.QJfor the 
considered i?-range. The fitted parameters are, for each 
sample, the Zr-Ol distance R\, the associated DW factor 
CT1 2 and AEo. The number of 01 neighbors, N±, is fixed 
to 6. The analysis of the EXAFS data for BaZr0 3 at 
11 and 300 K allows to determine a So 2 value equal to 
1 ± 0.07. This parameter was then fixed to 1 for all the 
further analysis. The EXAFS oscillations were refined in 
the i?-range [ 1.14-2.33 A]: the data for R lower than 1 A 
are affected by the background subtraction procedure^! 
and had to be excluded from the fit. Fits of good quality 
could be obtained, with reliability factors less than 1%, 
for all samples at all temperatures (Fig. [3] and Table [I]). 

In BaTio.65Zro.35O3 and BaZr0 3 samples, no signif- 
icant temperature dependence of the Zr-Ol distance is 
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FIG. 3: Modulus and imaginary part of the FT of k x(k) 
at 300 K. The dots represent measured data and the solid 
lines their best fit, obtained with the parameters values given 
in Table H] Measured data for R < 1 A are affected by the 
background subtraction procedure and were excluded from 
the analysis. 



TABLE I: Structural parameters deduced from the EXAFS 
analysis in the i?-range [ 1.14-2.33 A], for BaZr0 3 and BTZ 
relaxors. Ri is the length of the bond between an absorbing 
Zr atom and its first oxygen neighbor, and ai 2 the associated 
Debye- Waller factor. RF is the reliability factor of the fit. 
The uncertainties are of the order of ± 0.01 A and ± 0.0004 
A 2 for Ri and <ti 2 respectively. The energy shift AEq was 
found equal to 1.4 ± 0.8 eV for all fits. 
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detected, within the experimental accuracy of the EX- 
AFS technique. More surprisingly, the Zr-Ol distance 
hardly varies with the Zr substitution rate x, and keeps 
values very close to that found in BaZrC>3. 

The refined values of the a\ 2 DW factors in BaZrOa at 
11 and 300 K are in good agreement with those reported 
in Ref. 0. In BTZ relaxors, systematically higher a± 2 
values are obtained. From the definition of the DW fac- 
tor, this increase depicts either a stronger static disorder, 
or enhanced vibration amplitudes (associated with a de- 
crease of the Zr-Ol bond stiffness). As a matter of fact, 
the static and dynamic contributions can be separated, 
by analyzing the thermal evolution of the measured DW 
factor. In the absence of a static disorder, the DW factor 
only accounts for thermal vibrations and can be written, 
in the Einstein model, as: 

(T) = 2^iW 0th (i>' (2) 

where kg is the Boltzmann constant, and Mr the re- 
duced mass of all atoms involved in the scattering path. 
The Einstein temperature &e, characterizing the bond 
strength, increases with the bond stiffness. Given our 
hypothesis of a perfect Pm3m perovskite local structure 
in BaZrOs, the two <7i 2 values measured at 11 and 300 
K in this sample were fitted to Eq. [2J An Einstein tem- 
perature of 606 K ± 14 K was then obtained for the 
Zr-Ol bond. In BaTio.65Zro.35O3, the thermal evolu- 
tion of <j\ 2 is consistent with the same Oe value, pro- 
vided that a constant Aci 2 = + 0.0007 A 2 is added 
to Eq. [21 (Fig. This result gives evidence for the 

existence of a static distribution of the Zr-Ol distances 
in BaTio.65Zro.35O3, which is temperature-independent 
within the experimental accuracy. The full width at half 
maximum of a gaussian distribution of distances would 
be equal to 2\/21n2Aeri 2 = 0.06 A. This small value jus- 
tifies a posteriori the use of only one SS path in Eq. ^ 
the DW factor <7i 2 taking into account the small disorder 
that cannot be resolved in i?-space. 

The mean square deviation of the Ri distance linked 
to a static disorder, Acr^ 2 , is intended to describe a 
gaussian static distribution of distances in Eq. ^ with 
a 2 — o- 2 therm (T) + Aai 2 . However, a measured non-zero 
value of Aa 2 can describe other distance distributions 
(discrete or continuous), provided their width remains 
small. In the present case, the distribution of the Zr-Ol 
distances could be due to a distortion of the ZrOe units, 
coming from the Zr/Ti chemical disorder, and/or to a 
displacement of Zr atoms in their octahedral cages. To 
evaluate the magnitude of such a displacement, we cal- 
culated the FT of the EXAFS signals for Zr atoms in 
a perfect octahedron, weakly displaced either along the 
[100], [110], or [111] cubic axis. We then performed fits on 
the calculated signals in the same conditions as described 
above for the measured signals. From these simulations, 
we conclude that the refined DW factors presented in 
Table [I] for BaTio.65Zro.35O3 could correspond to a Zr 
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FIG. 4: Thermal evolution of the measured DW factor for 
the Zr-01 bond (symbols). The solid line represents Eq. |5] 
with 6e = 606 K. The dashed line corresponds to the same 
function shifted by 0.0007 A 2 . 



TABLE II: SS and MS paths used to analyze BaZrOs EXAFS 
data. Ni is the degeneracy of path i, and Ri its half-length, 
expressed as a function of Ri assuming an ideal perovskite 
structure. Zr c is the central absorbing Zr atom. The Ol 
and Ol' atoms are located on opposite sides of the Zr c atom, 
between Zr c and one of its six Zr third neighbors. 
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displacement of 0.07 A, in any of the three directions 
cited above. 

A static distribution of Zr-Ol distances is also expected 
in the two other relaxor samples, BaTio.70Zro.30O3 and 
BaTio.75Zro.25O3, from the txi 2 values measured at 300 K 
and by analogy with the case of BaTio.65Zro.35O3. 

In conclusion, the length and strength of the bond 
between a Zr atom and its first oxygen neighbor are 
found to be independent of the Zr substitution rate 
in BaTii-^Zr^Os samples. Only a small temperature- 
independent distribution of the Zr-Ol distances is de- 
tected in BaTio.65Zro.35O3, in opposition to the single 
distance in the regular octahedron present in BaZr03. 
These fluctuations of the Zr-Ol distance can be due to 
Zr displacements in their octahedra, but can also result 
from distortions of the ZrOe octahedra induced by the 
Zr/Ti chemical disorder. 

C. Analysis of the further neighbor contribution: 
BaZr0 3 

For R values up to 4.5 A, the FT of the EXAFS sig- 
nal involves scattering processes beyond the first oxygen 
neighbors. The amplitudes and phases of all the possi- 
ble paths have been calculated by using the FEFF8.00 
code, for an atomic cluster of 145 atoms representative 
of the BaZrOs Pm3m perovskite structure. The scatter- 
ing processes with relative weight lower than 2.5 % have 
been neglected. The single backscattering processes on 
the second (Ba), third (Zr), and fourth (02) neighbors 
of the Zr central atom must be considered, as well as 
several collinear MS paths and only two non-linear MS 
paths. One of the two latter paths involves Ba atoms. 
The other one is a triple scattering path within the ZrOe 
octahedron, which corresponds to a i?-range between the 
two main peaks of the FTs. We verified that it does 



not affect the fit, and did not consider it further. The 
retained paths for the EXAFS analysis are listed in Ta- 
ble [D] Note that the half-lengths of all the collinear MS 
paths are equal to the Zr-Zr distance. As a consequence, 
these MS paths contribute to the FT in the same R- 
range as the single backscattering path by a Zr atom, 
which makes the analysis complex. 

Assuming a perfect Pm3m cubic structure, all Ni are 
known and the Ri parameters for the ten paths can be 
expressed as a function of i?i , the distance between the 
central absorbing Zr atom and Ol (see Table ITT|) . The 
number of ai 1 parameters needed could be decreased us- 
ing the expression of the DW factor 38 : 

CT * 2 = 4^E^' +1 ~ ^)-^jj"+i] 2 )«' ( 3 ) 

3 

where t is the time, j represents an atomic site of the 
scattering path i, itj is the displacement vector of the j 
atom, and Kjj+i is the directing unit vector between j 
and j + 1 atoms at equilibrium. The resulting relations 
between the ai 1 parameters in BaZrOs are given in Table 
UTI For the fit of the measured FT, in the i?-range [ 1.14- 
4.52 A], seven parameters were then refined: AEq, 
C2 2 , 03 2 , (T4 2 , eg 2 , and trio 2 - Although they have already 
been determined in Sec. IIII Bl R\ and the strongly corre- 
lated parameter AEq were refined^, all distances up to 
4.5 A being functions of R\. On the other hand, the Sq 2 
and <ti 2 parameters were fixed to the previously refined 
values. The obtained parameters at 11 K and 300 K are 
presented in Table IIIII During the refinement process, 
the (T4 and erg 2 values were found to have no influence 
on the RF factor. These two DW factors are related to 
paths 4, 6, and 9, which give rise to very large, and rela- 
tively weak contributions in i?-space. Consequently, 04 2 
and erg 2 values do not affect significantly the amplitude 
and shape of these contributions and cannot be precisely 
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FIG. 5: Modulus and imaginary part of the FT of k 2 \{k) 
for BaZr0 3 , at 11 K (top) and 300 K (bottom). The dots 
represent measured data and the solid lines their best fit. The 
fit parameters are given in Table ITTTl 



TABLE III: Structural parameters deduced from the EXAFS 
analysis in the ii-range [ 1.14-4.52 A], for BaZrOg. The pa- 
rameters have already been defined in Tables HI and ITT1 The 
uncertainties are of the order of ± 0.01 A for Hi, ± 0.0003 
A 2 for (7 2 2 and <t 3 2 , and ± 0.003 A 2 for <7io . AEq was found 
equal to 0.0 ± 0.5 eV for both fits. 

Ri cr2 2 o- 3 2 crio 2 RF 

(A) (A 2 ) (A 2 ) (A 2 ) factor (%) 

300 K 2.10 0.0072 0.0048 0.013 1.04 

11 K 2.10 0.0026 0.0025 0.006 1.01 



determined, so that they are not given in Table ITTTl 

The calculated FTs are in very good agreement with 
the measured ones (Fig. EJ). The DW factors for the Zr- 
Ba (a-2 2 ), Zr-Zr (CT3 2 ) and Zr-02 (gin 2 ) bonds are con- 
sistent with those reported in Ref. |43[ Their thermal 
evolution was fitted to Eq. [21 resulting in Einstein tem- 
peratures equal to 194 ± 8 K, 263 ± 10 K, and 296 ± 50 
K for the Zr-Ba, Zr-Zr, and Zr-02 bonds respectively. 
The analysis of the EXAFS oscillations in BaZrC>3 
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FIG. 6: Modulus of the FT of the individual contributions 
to k 2 x(k), in BaZr0 3 at 300 K. SS paths, MS paths within 
Zr c octahedra, and MS paths outside Zr c octahedra are rep- 
resented as solid, dash-dotted, and dotted lines respectively. 
Numbers correspond to the path indexes of Table ITTTl 



gives information on the relative contributions of the dif- 
ferent paths to the FTs in the f?-range [ 2.5-4.5 A], which 
will help to analyze BTZ samples (see Fig. EJ. Only the 
backscattering path on Ba atoms (path 2 in Table ITT|) 
contributes in the i?-range [3.0-3.6 A], which allows a 
precise determination of the Zr-Ba distance and the as- 
sociated DW factor. The contributions in the i?-range 
[3.6-4.2 A] mainly arise from the MS linear paths that 
involve the third neighbor Zr atoms (paths 5 and 8 in 
Table [nj . The backscattering path on Zr atoms (path 
3) also contributes in the same i?-range, but it is much 
less important. The remaining MS paths (paths 4, 6, 7, 
and 9 in Table [H| correspond to weak contributions to 
the FTs, but it is necessary to take them into account 
in order to obtain good quality fits. The DW factors of 
paths 4, 6, and 9 were shown to have no influence on the 
fit. Finally, the backscattering path on 02 atoms (path 
10) contributes over a very large i?-range [3-4.6 A]. It is 
then compulsory to keep this path in the analysis, as it 
cannot be separated from the other contributions to the 
FTs. 



D. Analysis of the further neighbor contribution: 

BaTii„ a; Zr :E 03 relaxors 

In the following analysis, we shall introduce two differ- 
ent models of the local structure in BTZ relaxors. In a 
first step, we show that a model with aligned Zr, 01, and 
(Zr/Ti) atoms in the perovskite structure cannot account 
for the measured EXAFS signals. In a second step, we 
introduce a buckling angle for the Zr-01-(Zr/Ti) bonds. 
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1. Basic model for the BTZ analysis 



The BaZr03 model (see Sec. IIII CI) is used as a start- 
ing point to analyze BTZ samples. To take into account 
the Ti third neighbors, we also calculated the amplitudes 
and phases of the scattering paths within a hypothetical 
cluster of BaTiOs in the Pm3m cubic perovskite struc- 
ture, with a single Zr impurity as the absorbing atom. 
For the paths that do not involve Ti atoms, the calcu- 
lated amplitude and phase remain strictly identical to 
those calculated for the BaZrC-3 structure. Thus, we 
only consider three extra paths 3' (Zr c — * Ti — ► Zr c ), 
5' (Zr c -► Ti -> 01 -> Zr c ), and 8' (Zr c -> 01 -► Ti -> 
01 — * Zr c ), in addition to those presented in Table ITT1 
The fc-dependence of these path amplitudes differs sig- 
nificantly from that of the corresponding paths 3, 5, and 
8 of the BaZr0 3 model. Thus, the EXAFS technique al- 
lows the refinement of structural parameters for Zr and 
Ti atoms with limited correlations effects, despite similar 
Zr-Zr and Zr-Ti distances. A new parameter Nz r has to 
be introduced, which is the average number of Zr third 
neighbors of the central atom Zr c . The degeneracy of 
paths 3, 5, and 8 is then multiplied by Nz r /6, and that 
of paths 3', 5', and 8' by (6-iVz r )/6. A unique ay 2 pa- 
rameter is attached to the paths 3', 5' and 8', according 
to Eq. 

Zr atoms being larger than Ti ones, they likely cause 
local distortions in the BaTi03 matrix. This difference in 
size is taken into account by introducing independent Zr- 
01, Zr-Ba, Zr-Ti, and Zr-02 distances. The 01 atoms ly- 
ing between Zr c and (Zr/Ti) atoms in this first model, the 
Zr-Zr and Ti-01 distances are defined using the expres- 
sions dzr-Zr = 2dzr-01 and dri-Ol = dzr-Ti — dzr-Ol- 

The DW factors 04 2 and erg 2 were fixed to the same values 
as in BaZr03, since they are not expected to affect the 
fit (see Sec. IIII C|) . Furthermore, the relations between 
the DW factors presented in table ITT1 remain correct, so 
that only er 2 2 , cr 3 2 , er 3 ' 2 , and cio 2 DW factors were re- 
fined. The Zr-01 distance (i?i) and ci 2 , as well as AEq 
and So 2 , were fixed to the values refined in the i?-range 
[1.14-2.33 A] (see Table0. 

Within this model, the fitting procedure over the R- 
range [1.14-4.52 A] yields non-physical values of the CT3 2 
and ay 2 parameters (0.033 and -0.0008 respectively, at 
11 K). In fact, the calculated contribution of paths linked 
to Zr atoms (paths 3, 5, and 8 in Table [HJ are com- 
pletely damped (see top of Fig. 01. Only paths involv- 
ing Ti atoms then contribute to the calculated signal, 
in the i?-range [3.4-3.9 A]. Fixing 03 2 to its value in 
BaZrOs made impossible a good agreement between cal- 
culated and experimental signals. In particular, the shift 
of the FT's imaginary part in the i?-range [3.6-3.9 A] in- 
dicates an overestimated Zr-Zr distance (see bottom of 
Fig. UJ. Hence, wc deduce that dzr-Zr must be smaller 
than 2dzr-oi- 
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FIG. 7: Modulus and imaginary part of the FTs of 
the observed (dots) and refined (solid line) k 2 x(k) for 
BaTio.65Zro.35O3 at 11 K. The modulus of the FT for paths 
3 (dashed line) and 3' (dotted line), which describe the 
backscattering on Zr and Ti atoms respectively, are also rep- 
resented. The fits were obtained in the hypothesis of aligned 
Zr c , Ol and (Zr/Ti) atoms. Top: the path 3 is completely 
damped by a non-physical, huge 03 2 DW factor. Bottom: 03 2 
is fixed to its value in BaZrC>3, i.e. 0.0020 A 2 . The calculated 
imaginary part is then shifted towards high _R-values in the 
ii-range [3.7-3.9 A], where the path 3 now contributes. The 
latter clearly indicates an overestimated dzr-Zr distance. 



2. Determination of the buckled local structure in BTZ 
relaxors 

The hypothesis of aligned Zr c , Ol and Zr atoms be- 
ing not valid in BTZ samples, a new model has to be 
built, where the Ol atoms are no longer aligned with Zr 
ones. In order to get all different possible configurations, 
one can introduce two mean buckling angles, Qzr and 
9 Ti , defined as 180° - ZrOlZr and 180° - ZrOlTi re- 
spectively. All the possible paths and their parameters 
were calculated in the preceding BaZr03 and BaTi03 
clusters, moving the Ol atoms out of their collinear sites 
with varying values of the buckling angle O. The only 
effect of the Ol atom displacements is the damping of the 
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FIG. 8: Evolution of the focusing path amplitudes with in- 
creasing values of the buckling angle O: example of path 5. 



amplitudes of the MS paths in which these atoms play a 
focusing role, i.e. paths 5, 5', 8, and 8' ( Table ITTf) . As an 
example, we show the evolution of the amplitude of path 
5 as a function of on Fig. [S] This evolution can be 
expressed by expending the scattering amplitude of the 
focusing paths about 9 = 0: 

A(fc,e)w A(fc,o)[i-&(fc)eT, (4) 

where n is the number of scattering processes by the off- 
centered focusing atom (01)^. For each of the paths 5, 
5', 8, and 8', the coefficients b(k) were determined for 
various buckling angles and averaged. Eq. Q]is found to 
be correct for buckling angles values below 20°. 

The expressions for Zr-Zr and Ti-Ol distances are 
now given by: dz r -Zr — 2dz r _oi C0S (©Zr/2), and 

dn-oi = \J d 2 Zr _ Ti - d 2 Zr _ 01 sin 2 Q Ti - d Zr -o\ cos <d Tl , 
from which the lengths of paths 1 to 8 can be easily de- 
rived. The path length Rg splits into several different 
values, due to the misalignment of the 01 oxygen atoms. 
This effect can be modeled by an increase of the cor- 
responding DW factor erg 2 , but we already showed that 
this parameter does not affect the fit quality. Thus, we 
determined the mean distance Rg with the expression 
given in Table [H] using the the erg 2 value of BaZrC^. 
The length of path 10 splits as well, depending on the 
size and shape of the octahedra next to the central one 
Zr c Og- This information being out of the scope of this 
study, we assigned to this path an average value of the 
Zr-02 distances, i?io- In that case, 0io 2 is expected to 
take into account both thermal and static disorders. In 
the absence of an adequate model for the 02 reparti- 
tion, the values obtained for i? 10 and a w 2 will not be 
interpreted. Concerning the DW factors, the linearity of 
paths 5, 5', 8 and 8' is broken by the off-centering of 01 



atoms with respect to Zr-(Zr/Ti) bonds, so that the re- 
lations (75 2 = (7 8 2 = (73 2 and 05 = Og = cr 3 / 2 are no 
longer exact. However, we postulated their validity for 
the buckled paths, in order to limit the number of free 
parameters. The buckling angles <d Zr an d Q?i were de- 
termined separately, as preliminary fits using the same 
buckling angle for both Zr c -01-Zr and Zr c -01-Ti bonds 
did not yield a good agreement between experimental and 
calculated signals (see Fig. |5} ■ Ten parameters were then 
refined: R 2 , i? 3 ', -Rio, N Zr , Q Zr , Ti , er 2 2 , er 3 2 , a 3 , 2 , and 
erio 2 . The (73 2 and ay 2 parameters were found to reach 
very low (non- physical) values when not fixed. Like N Zri 
0Ti ; an d ®zr, these two DW factors determine the signal 
amplitude in the i?-range [ 3.4-3.9 A]. Strong correlation 
effects exist between these five parameters and various 
sets of their values give rise to similar, very good fits. 
To go further in the analysis, we assume that <73 2 and 
cr-y 2 in BTZ relaxors are equal to cr 3 2 in BaZrOs, which 
is determined for all temperatures using Eq. [21 with the 
Einstein temperature of 263 K reported in Section IlII CI 
Values of refined parameters in this hypothesis are pre- 
sented in Table IIVI and fits are shown on figure ^3 f° r 
the three BTZ relaxors. The fits are not very sensitive 
to the variations of ©Ti in the range [0-10°]. ©Ti values 
larger than 10° are excluded since they do not yield sat- 
isfactory agreements. Note that the Zr-Ba distance and 
the associated DW factor 02 2 are accurately defined, the 
contribution of path 2 to the FTs being well separated 
from the others. 

One needs to estimate the consequences of fixing 03 2 
and 03' 2 for the analysis of BTZ relaxors. For this pur- 
pose, we performed several fits for various, not necessarily 
equal, values of 03 2 and 03' 2 , in a range of ± 0.0010 A 2 
around the CT3 2 value in BaZrOs. The Q Zr , N Zr , and Ry 
values were found to vary by ± 1.5°, ± 0.2, and ± 0.01 A 
respectively. Furthermore, the refined values of & Zr be- 
ing at the limit of the validity range of Eq. 0] one expects 
an error of ± 1° on its determination. The uncertainties 
given in Table El are those determined by the FEFFIT 
code, augmented by the preceding amounts. Note that in 
the case of BaTio.65Zro.35O3, the refined value of N Zr is 
found to be independent of the temperature, as it should 
be. This confirms a posteriori that the hypothesis of an 
identical thermal evolution of 03 2 and 03' 2 in BTZ sam- 
ples and 03 2 in BaZrOs is reasonable. 



3. Fit results 

For all the relaxor samples investigated, the refined 
values of the buckling angle <d Zr range from 18 to 20°. 
©Ti being smaller (less than 10°), the ZrOg octahedra 
are therefore distorted differently depending on the type 
of the third neighbors (Zr or Ti) of the Zr central atom. 

In BTZ relaxors, the existence of a nearly constant 
buckling angle whatever the Zr substitution rate x, to- 
gether with the invariant Zr-Ol distance (Sec. IIII Bp re- 
sult in a Zr-Zr distance of the order of 4.15 A, which 
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TABLE IV: Structural parameters deduced from the EXAFS analysis in the _R-range [ 1.14-4.52 A], for BTZ relaxors. The 
path half-lengths are denoted Ri, and the associated DW factors a 2 . Qz r is the buckling angle of the Zr-Ol-Zr bonds, defined 
as 180° - ZrOlZr, and Nz r the mean number of Zr third neighbors of the Zr c atom. 7?3 is deduced from R\ (Table^J and Qzr- 
The uncertainties are of the order of ± 0.006 A, ± 0.0004 A 2 , ± 5°, ± 0.04 A, ± 0.02 A, ± 0.04 A, ± 0.005 A 2 , and ± 1.0 for 
R2, cr 2 2 , Qzr, R3, Rs', Rio, o"io 2 and N Zr respectively. 
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FIG. 9: Modulus and imaginary part of the FTs of 
the observed (dots) and refined (solid line) k 2 x(k) for 
BaTio.65Zro.35O3 at 11 K. The modulus of the FT for paths 
3 (dashed line) and 3' (dotted line), which describe the 
backscattering on Zr and Ti atoms respectively, are also rep- 
resented. We show the best fit obtained in the hypothesis of 
a single buckling angle for Zr c -01-Zr and Zr c -01-Ti bonds, 
for reasonable values of 03 2 and a 3 2 parameters. The miss- 
ing amplitude in the ii-range [3.3-3.7 A], characteristic of the 
path 3', indicates an overestimated buckling angle for the Ti 
bonds. 



hardly varies with x. Of course, the Zr-Zr distance is 
smaller in BTZ samples than in BaZrOa, since the oxy- 
gen atoms are aligned with the Zr-Zr bonds in the latter 
compound. The Zr-Ti distance (4.07 A) does not change 
with x, and is found to be significantly smaller than the 
Zr-Zr distance. The Zr-Ba distance, on its turn, increases 
with x. From the Zr-Ti distance and the buckling angle 
of the Zr-Ol-Ti bond, we can derive the Ti-Ol distance 
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FIG. 10: Modulus and imaginary part of the FT of k 2 \{k) at 
either 11 K or 300 K, for BTZ samples. The dots represent 
measured data and the solid lines their best fit, obtained with 
the parameters values given in Table HVI 



value at 300 K in BTZ relaxors, which is of the order of 
1.98 A. This corresponds to the shortest Ti-0 distance 
in the quadratic BaTi0 3 at 300 K^. 

We now consider the temperature dependence of 02 2 in 
BaTio.65Zro.35O3 sample. As shown on Fig. 1111 it can be 
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FIG. 11: Thermal evolution of the measured DW factor for 
the Zr-Ba bond (symbols). The solid line represents Eq. [5] 
with 6e = 194 K. The dashed line corresponds to the same 
equation with 6 E = 231 K, and shifted by 0.0019 A 2 . 



described following Eq. [2] with an Einstein temperature 
equal to 231 ± 4 K, instead of 194 ± 8 K in BaZr0 3 . 
This result indicates an increase of Zr-Ba bond strength 
with increasing Ti-content. In addition, a constant A02 2 
= 0.0019 A 2 must be added to Eq. [2J revealing the pres- 
ence of static disorder. The full width at half maximum 
of the corresponding gaussian distance distribution would 
be equal to 2\/2 In 2A02 2 = 0.10 A. However, let us re- 
mind that other distance distributions can be expected 
(see Sec. MI Bl) . The possible distribution of the Zr-Zr 
and Zr-Ti bond lengths is out of the scope of this study, 
since 03 2 and a^ 2 parameters had to be fixed during the 
refinement process. 

The mean number of Zr third neighbors around a Zr 
atom Nzr brings information on the distribution of the 
Zr and Ti atoms in BTZ relaxor samples. Nz r is found 
to be equal to 2.7, 3.1, and 3.4 in BaTio.75Zro.25O3, 
BaTio.70Zro.30O3, and BaTio.65Zro.35O3 samples respec- 
tively. Despite the large uncertainty attached to Nz r 
(± 1.0), the measured values are higher than those ex- 
pected for a random distribution of Zr and Ti atoms: 
1.5, 1.8, and 2.1. This observation indicates a tendency 
of Zr atoms to segregate in Zr-rich regions. However, 
the EXAFS technique giving only an average of the local 
structure, we cannot conclude on the Zr concentration 
in the Zr-rich regions, and hence on the sizes of these 
regions. 

In summary, in all the BTZ relaxors investigated the 
shape of a ZrO@ octahedra depends on the nature of 
its neighboring octahedra (ZrOg or TiOg), the oxygen 
atoms being differently off-centered with respect to the 
Zr-Zr and Zr-Ti bonds. The Zr-Zr and Zr-Ti distances 
do not depend on x, Zr-Ti distances being significantly 
shorter than Zr-Zr ones. The Zr-Ba bond length increases 



with the Zr-content, and is affected by a strong static 
disorder. Finally, the EXAFS analysis reveals a ten- 
dency of Zr atoms to segregate. Note that in the case of 
BaTio.65Zro.35O3, the Zr environment is found to be re- 
markably stable with temperature. Therefore, no impor- 
tant changes of the local structure are observed around 
Zr atoms in the temperature range of the maximum of 
the dielectric permittivity. 



IV. CONCLUDING DISCUSSION 
A. Microstructural picture of BTZ relaxors 

In the preceding sections we have described quantita- 
tively the local structural features of the Zr environment 
in BTZ relaxors. It is now interesting to discuss how this 
local structure contrasts with the average, long-range cu- 
bic structure evidenced by X-ray diffraction (XRD) . The 
interatomic distances deduced from EXAFS are com- 
pared to those determined from XRD on Fig. ^] The 
two techniques yield quite different values of distances. 
On one hand, the XRD measurements show that the unit 
cell dimensions of BTZ samples follow the Vegard's law, 
i.e., the cell volume linearly increases from its value in 
BaTi0 3 (64.286 A 3 ) to that in BaZr0 3 (73.665 A 3 ). In 
the cubic relaxor samples, all the average distances then 
linearly increase with the Zr substitution rate x. On the 
other hand, the distances deduced from the EXAFS anal- 
ysis systematically exceed those expected in the average 
structure and, with the exception of the Zr-Ba distance, 
are found to be independent of x. Considering that EX- 
AFS probes the structure on a very local scale, the latter 
observation is direct evidence that the local structure is 
different from the average structure. 

We have shown in Sec. MI D 31 that Zr-atoms tend 
to segregate in BTZ relaxors. The latter result is in 
agreement with previous speculations based on Raman 
scattering 23 . If we assume that the Zr-rich regions con- 
sist of BaZr03 spherical inclusions in a BaTi03 matrix, 
one can calculate the mean number of segregated atoms. 
It is estimated to 27, 15, and 10 in BaTio.65Zro.35O3, 
BaTio.70Zro.30O3, and BaTio.75Zro.25O3 respectively, 
which corresponds to diameters of 14.1, 11.5, and 10.8 A. 

The EXAFS study of BTZ relaxors at the Zr K-edge 
allowed us to evidence local structural deviations from 
the average cubic structure, as well as chemical inhomo- 
geneities. In the following, we shall propose a microstruc- 
tural picture of BTZ relaxors and discuss its origin and 
implications on the relaxor behaviour. In order to ex- 
pose the possible basic physical mechanisms in a clear 
and simple manner, we will assume that the segregation 
of Zr atoms consists of small inclusions of BaZr03 in a 
BaTi03 matrix. However, the real repartition of the Zr 
atoms could be much more complex (see Sec. MI D 3ll . 

Such BaZr03 inclusions in BTZ relaxors must be sub- 
mitted to a chemical pressure from the BaTi03 environ- 
ment, which has a smaller unit cell volume than BaZr03. 
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polyhedra more rigid than AO12 ones 51 . Furthermore, 
it has been proposed in the literature that BaZr03 un- 
dergoes a pressure-induced phase transition to the tilted 
a~a~a~ perovskite structural (Glazer's notations^). 
Assuming that the BaZi'03 inclusions adopt the latter 
structure in BTZ relaxors, the tilt angle calculated from 
the buckling angle of the Zr-O-Zr bonds is about 11°, 
which is a reasonable value for a tilted perovskite^S*^. 
In summary, within our simplified picture, BaZrC>3 inclu- 
sions do not adopt the cubic bulk structure of BaZrC>3 
but present considerable distortions away from it. The Zr 
environment in BTZ relaxors appears to be determined 
by both the high stiffness of the ZrOe octahedra and 
chemical pressure effects. 

Let us now discuss the structure around the inclusions. 
Given that the Zr-Zr distance is approximately 2% larger 
than the average B-B cation distance deduced from XRD 
(Fig. 1121-c). we expect a significant strain at the interface 
of the BaZrC>3 inclusions and the BaTi03 matrix. The 
TiC>6 octahedra adjacent to the inclusions likely accom- 
modate the structural misfit between BaZrC>3 inclusions 
and the BaTiC>3 matrix. Considering that the cell vol- 
ume of BaZrC>3 is larger than that of BaTiC>3, we expect 
that these adjacent TiC>6 octahedra are submitted to a 
tensile strain. Although their deformation shape remains 
an open question, different Ti + displacements (and thus 
polarity) are expected in the adjacent regions when com- 
pared to the BaTiOs matrix. 



FIG. 12: Evolution of the Zr-Ol (a), Zr-Ba (b) and Zr-(Zr/Ti) 
(c) distances with the Zr substitution rate x. EXAFS results 
are reported as round dots. They are compared to the cor- 
responding average distances deduced from the X-ray diffrac- 
tion experiments (crosses); i.e. a/2, aV3/2, and a respec- 
tively, where a is the cubic cell parameter. 



It is first interesting to note that the Zr-0 distance in 
BTZ relaxors is almost insensitive to this chemical pres- 
sure, since it is the same as in BaZrOs (Fig. 1121a). On 
the other hand, the Zr-Ba distance decreases with in- 
creasing Ti-content (Fig. 1121-h) and hence, is more sen- 
sitive to the chemical pressure. This latter difference in 
ABO3 perovskites is well understood within the polyhe- 
dral approach: the A-0 distance in perovskites is gen- 
erally more compressible than the shorter B-0 distance 
in the closed-packed BO@ octahedra™. Insight into the 
BOq linkage can be also gained from the ZrOZr angle 
value, which deviates considerably from 180° C Table llVjl . 
In a perovskite structure, the presence of buckled Zr-O- 
Zr bonds suggests the presence of octahedron tilts. This 
distortion could be a direct consequence of the chemi- 
cal pressure exerted by the surrounding BaTi03 matrix. 
Within this context, let us remind that hydrostatic pres- 
sure reveals and/or enhances structural lattice instabili- 
ties related to tilts in the perovskite compounds with -806 



B. Comparison to PZT 

We have already mentioned in the introduction that 
the ferroelectric PbZri^Ti^Os (PZT), similar to BTZ, 
does surprisingly present no relaxation, making the com- 
parison of BTZ to PZT insightful. First, we note that in 
both BTZ and PZT, the mean Zr-0 distance is almost 
unaffected by the Ti-content, and the observed Zr-0 dis- 
tance distribution is of the same order of magnitude^. 
As a consequence, the relaxor properties of BTZ cannot 
be linked to the first neighbour environment of Zr atoms. 

A key feature for understanding the different behaviors 
of BTZ and PZT could lie in the local chemical order, 
which is known to play a determinant role for relaxor 
properties (e.g. the so-called 1:1 Mg/Nb order on the 
nanometer scale in PMN)i2iiI. In that context, it is in- 
teresting to note that BTZ cannot be synthesized in the 
region 0.5 < x < 1 (Ref. 0), which suggests an important 
internal chemical strain which inhibits the formation of 
BTZ across the whole composition range. We speculate 
that the local phase segregation observed for BTZ relax- 
ors is a precursor for the macroscopic phase separation 
observed for higher Zr substitution rates. Such an in- 
ternal chemical strain does not seem to exist for PZT 
since it forms a solid solution whatever the Zr/Ti ra- 
tio. The distribution of elastic fields is thus expected to 
differ significantly from BTZ to PZT. In particular, the 
strain around the local phase separation in BTZ can be 
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regarded as the source of the elastic random fields which 
have been proposed earlier by Farhi et al2i. Such elastic 
random fields (similarly to electric random fields in com- 
mon relaxors) are then considered to lead to particular 
local pattern of polar regions, which become the source 
of relaxor properties. 
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